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We present studies of the time-evolution of a two- 
component system of Bose-Einstein condensates (BEC) in the 
\F = l,mf = —1} and |2, 1} spin states of 87 Rb. The two con- 
densates are created with a well-defined relative phase and 
complete spatial overlap. In subsequent evolution they un- 
dergo complex relative motions that tend to preserve the to- 
tal density profile. The motions quickly damp out, leaving 
the condensates in a steady state with a nonnegligible (and 
adjustable) overlap region. 



Since its realization in dilute atomic gases , Bose- 
Einstein condensation (BEC) has afforded an intriguing 
glimpse into the macroscopic quantum world. Attention 
has recently broadened to include exploration of systems 
of two or more condensates, as realized in a magnetic trap 
in rubidium [Q and subsequently in an optical trap in 
sodium H . Theoretical treatment of such systems began 
in the context of superfluid helium mixtures || and spin- 
polarized hydrogen [7], and has now been extended to 
BEC in the alkalis g^. 

The first experiments involving the interactions be- 
tween multiple-species BEC were performed with atoms 
evaporatively cooled in the \F = 2, to/ = 2) and |1, —1) 
spin states of 87 Rb fl. These experiments demonstrated 
the possibility of producing long-lived multiple conden- 
sate systems, and that the condensate wavefunction is 
dramatically affected by the presence of interspecies in- 
teractions. In this Letter, we report results from initial 
studies of simultaneously trapped BECs in the |2, 1) and 
|1, — 1) states of 87 Rb (denoted hereafter as |2) and |1), 
respectively). The two states are completely distinguish- 
able since the hyperfine splitting is much larger than any 
other relevant energy scale in the system. We produce 
arbitrary superpositions of |1) and |2) that begin with a 
well-defined relative phase, spatial extent, and "sag" - 
the position at which the magnetic trapping forces bal- 
ance gravity for each state. The fine experimental control 
of this double condensate system permits us to study its 
subsequent time-evolution under a variety of interesting 
conditions, most notably those in which there remains 
substantial overlap between the two states. 

The apparatus and general procedure we use to at- 
tain BEC in Rb are identical to those of our previous 
work ||ll| and will be reviewed here but briefly. We use a 



87 



rf photon 
(1-4 MHz) 



Rb 



■ |2,-2> 



|1,-1) 





|2> 



|2,2) 



F = 2 



microwave photon 
-6.8 GHz 



|1,0> 



F = 1 



|1,1> 



FIG. 1. A schematic of the ground-state hyperfine levels 
(F = 1, 2) of 87 Rb, shown with a splitting due to an applied 
magnetic field. The two-photon transition is driven between 
the 1 1,-1) (|1>) and |2, 1) (|2)) states. 



double magneto-optical trap system to load roughly 10 9 
|1) atoms into a time-averaged, orbiting potential (TOP) 
magnetic trap |jL2|j . The atoms are magnetically com- 
pressed and evaporatively cooled Q for 30 s until they 
form a condensate of approximately 5 x 10 5 atoms with 
no noticeable non-condensate fraction (we estimate that 
> 75% of the entire gas is in the condensate). After 
completion of the evaporation cycle, the magnetic trap 
is ramped adiabatically to various bias fields and spring 
constants for the subsequent experiments. 

The double condensate system is prepared from the 
single |1) condensate by driving a two-photon transi- 
tion pd| ] consisting of a microwave photon near 6.8 GHz 
and a radiofrequency (rf) photon of 1-4 MHz, depending 
on the Zeeman splitting (Fig. |I|). As in pdj ], we are able 
to transfer quickly any desired fraction of the atoms to 
the 1 2) state by selecting the length and amplitude of the 
two-photon pulse. The two condensates are created with 
identical density distributions, after which they evolve 
and redistribute themselves for some time T. We then 
turn off the magnetic trap and allow the atoms to expand 
for 22 ms for imaging. 

We selectively image the densities of either of the two 
states (ni and 712) or the combined density distribution 
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(ut) by changing the sequence of laser beams applied to 
the condensates for probing jll| . Since the expansion and 
imaging are destructive processes, each image is taken 
with a different condensate; the excellent reproducibility 
of the condensates permits us to study the time-evolution 
of the system by changing the time T. We subsequently 
use a least-squares algorithm to reconstruct the relative 
positions of the condensates from the images of n%, 112, 
and riT at each time T — a task made necessary by shot- 
to-shot jitter in the image positions on the CCD-array 
detector. 

The evolution of the double condensate system, in- 
cluding the release from the trap and subsequent ex- 
pansion |l],[l3|], is governed by a pair of coupled Gross- 
Pitaevskii equations for condensate amplitudes $i: 
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where i,j = 1,2 (i ^ j), Vi is the magnetic trapping 
potential for state i, the mean-field potentials are Ui = 
47r?i 2 a,;|$i| 2 /m and = 47r7i \$>j\ 2 /m, m is the mass 
of the Rb atom, and the intraspecies and interspecies 
scattering lengths are at and 0^. In the Thomas-Fermi 
limit, the condensate density distributions are dominated 
by the potential energy terms of Eq. (|]J) . Consequently, 
the expanded density distributions retain their spatial 
information and emerge with their gross features (such 
as the relative position of the condensates) intact. 

The similarity in scattering lengths Oi, a,2, and a%2 
implies that the total density ut will not change sig- 
nificantly from its initial configuration even though the 
two components may redistribute themselves dramati- 
cally during the evolution time T. In 87 Rb, the scattering 
lengths are known at the 1% level to be in the propor- 
tion ai : a\2 : a2 :: 1.03 : 1 : 0.97, with the average of 
the three being 55(3) A HQ]- The near-preservation of 
the total density tit can be approached theoretically by 
deriving from Eq. ([j]) the hydrodynamic equations of mo- 
tion [ fisl for tit and evaluating them in the limit that the 
fractional differences between the scattering lengths are 
small. The pressures that tend to redistribute ut must 
also be small. A similar argument pertains if the min- 
ima of the trapping potentials V% and V2 are displaced 
from one another (see below) by a distance that is small 
compared to the size of the total condensate; once again, 
the effects on the equilibrium distribution of the individ- 
ual components may be profound but the total density 
should remain largely unperturbed Jl6[ . 

The rotating magnetic field of the TOP trap gives 
rise to a subtle behavior that permits us to displace the 
minima of the trapping potentials V\ and V2 with respect 
to one another |]l7| , |l8| . In the rotating frame, the two 
states see two different magnetic fields as a function of 
the bias field rotation frequency and sense of rotation (as 
well as the strengths of the bias and quadrupole fields). 




FIG. 2. (Color) (a) The image of the jl) condensate ex- 
hibits a "crater," corresponding to a shell in which the |2) 
atoms (b) reside. For this trap, v z = 47 Hz with zero 
relative sag. By changing the strength of the magnetic 
quadrupole field we can introduce a nonzero relative sag, 
which shifts the location of the "crater" (c). (Each square 
in this post-expansion image is 136 fim on a side.) 



By adjusting these parameters, we can change the sign of 
the relative sag or cause it to vanish |L9| while preserving 
(to first order) the same radial (v r ) and axial [y z — y/8v r ) 
trap oscillation frequencies. 

In a first experiment, we choose a trap that has zero rel- 
ative sag (y z = 47 Hz) and transfer 50% of the atoms to 
the 1 2) state with a ~ 400 /is pulse. When T — 30 ms, we 
observe a "crater" in the image of the |1) atoms (Fig. |]a). 
The "crater" corresponds to a region occupied by the 
1 2) atoms (Fig. ||b), indicating that the |1) atoms have 
formed a shell about the |2) atoms. This is consistent 
with the theoretical observation that it is energetically 
favorable for the atoms with the larger scattering length 
(|1)) to form a lower-density shell about the atoms with 
the smaller scattering length (|2)) ||. At longer times 
the condensates separate from one another radially |^0| . 

In order to explore the boundary between the two con- 
densates, we perform a series of experiments in a trap in 
which we displace the trapping potentials such that the 
minimum of V2 is 0.4 (im lower than that of Vy , or approx- 
imately 3% of the (total) extent of the combined density 
distribution in the vertical direction. The subsequent 
time-evolution of the system is shown in Figs. |3 and |^. 
The two states almost completely separate (Fig.lSa-c) af- 
ter 10 ms; they then "bounce" back until at T = 25 ms 
the centers-of-mass are once more almost exactly su- 
perimposed (Fig. |]), although a distinctive (and repro- 
ducible) vertical structure has formed (Fig. [|d-e-f). By 
T = 65 ms, the system has apparently reached a steady 
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FIG. 3. (Color) Time-evolution of the double-condensate 
system with a relative sag of 0.4 /jm (3% of the width of 
the combined distribution prior to expansion) and a trap fre- 
quency v z = 59 Hz. 
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FIG. 4. The relative motion of the centers-of-mass of the 
two condensates under the same conditions as those in Fig. ^. 

state (Fig. §g-h, Fig. |J) in which the separation of the 
centers-of-mass is 20% of the extent of the cloud. From 
these images we observe: (i) the fractional steady-state 
separation of the expanded image is large compared to 
the fractional amount of applied symmetry breaking, as 
we expect for a repulsive interspecies potential; (ii) the 
placid total density profile (rightmost column of Fig. ||) 
betrays little hint of the underlying violent rearrange- 
ment of the component species; and (iii) the component 
separation is highly damped, although it is not yet cer- 
tain what mechanism is responsible. With respect 
to the damping, the excitation is in no sense small and 
may therefore be poorly modeled by theories that treat 
the low-lying, small-amplitude excitations |l0| of double 
condensates. 

Finally, we show the optical density as a function of 
relative number and position on the condensate vertical 
axis in order to better appreciate the amount of overlap 
between the two states at T = 65 ms (Fig. ||), which 
remains substantial despite the underlying separation. 
Each plot is averaged across a ~ 14 /im wide vertical 
cut through the centers of the two condensates. From 
the overlap shown, one could determine the magnitude 
of the interspecies scattering length a\2 by comparison to 
numerical solutions of the Gross-Pitaevskii equations (|l|) 
for our trapping conditions. Such a calculation is beyond 
the scope of the present work. 

In related work, we have read out the relative quantum 
phase of the condensates in the overlap region with a sec- 
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FIG. 5. Vertical cross-sections of the density profiles at 
T = 65 ms for different relative numbers of atoms in the 
two states. The combined density distribution (solid line) 
is shown for comparison to the Thomas-Fermi parabolic fit 
(dashed line). The trap parameters are the same as those in 
Fig. | 

ond two-photon pulse, making use of Ramsey's method 
of separated oscillatory fields. We expect the presence 
of a weak cw two-photon drive, resonant in the over- 
lap region, to result in phase-locking and phase-sensitive 
currents analogous to the Josephson effect in supercon- 
ductors [p2[ . We will discuss the phase-evolution of the 
condensates in a future publication. 
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